Entangled-like Chain Dynamics in Non-entangled Polymer Blends 
with Large Dynamic Asymmetry 
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We discuss simulations of a simple model for polymer blends in the framework of the Rouse model. 
At odds with standard predictions, large dynamic asymmetry between the two components induces 
strong non-exponentiality of the Rouse modes for the fast component. Despite chains being much 
shorter than the entanglement length, it also induces dynamic features resembling a crossover to 
entangled-like chain dynamics. This unusual behavior is associated to strong memory effects which 
break down the assumption of time uncorrelation of the external forces acting on the tagged chain. 

PACS numbers: 83.10.Rs, 83.10.Kn, 83.80.Tc 
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Polymer blends are systems of wide technological in- 
terest whose rheological properties can be tuned by 
varying the mixture composition. They are dynami- 
cally heterogeneous: starting from two homopolymers 
with different mobilities, two separated segmental re- 
laxations are still observed in the blend state. In the 
usual scenario the two components display qualitatively 
similar dynamic features in the blend. A rather differ- 
ent scenario arises, for low concentration of the fast 
component, if the two homopolymers exhibit very dif- 
ferent glass transition temperatures [3, • In this case 
the two components exhibit strong dynamic inmiscibil- 
ity, with a large separation (dynamic asymmetry) in 
their respective relaxation times. The latter can be of 
even 12 decades for high dilution of poly(ethylene ox- 
ide) (PEO) in poly (methyl methacrylate) (PMMA) Q. 
In such conditions the motion of the fast component 
seems to be strongly confined by the slowly relaxing ma- 
trix formed by the slow component. We have recently 
performed a computational investigation of the struc- 
tural a-relaxation in a simple bead-spring model for 
polymer blends ||] . The introduction of size disparity 
between the monomers of the two components induces 
a large time scale separation at low concentrations of 
the fast component, which displays unusual dynamic 
features, as logarithmic relaxation for density correla- 
tors, or dynamic decoupling between self-motions and 
intrachain collective relaxation. These unusual obser- 
vations are supported by recent fully atomistic simula- 
tions 0] and neutron scattering experiments 0] on the 
real blend PEO/PMMA. 

In this Letter we show that anomalous dynamic fea- 
tures induced by confinement on the fast component 
extend to large scales. By analyzing relaxation of the 
chain degrees of freedom within the framework of the 
Rouse model Q novel features are observed for the 
behavior of the chain normal modes of the fast com- 
ponent. We simulate short chains (4 < N < 21 



monomers) corresponding to the non-entangled case in 
the homopolymer state (the entanglement length [y] is 
A c ~ 35 monomers). This choice allows us to compare 
results between the blend and homopolymer states for 
chain lengths for which the homopolymer case shows 
reasonable agreement with predictions of the Rouse 
model @, 0, [1, @| ■ We test the two basic predictions: 
orthogonality and exponentiality of the Rouse modes. 
Simulations show that orthogonality is a reasonable ap- 
proximation for both components, as well as exponen- 
tiality in the case of the slow one. However, large devi- 
ations from exponentiality are found for the fast com- 
ponent in the blend. We observe a striking crossover, 
by increasing the dynamic asymmetry in the blend, to 
a regime resembling scaling features characteristic of 
entangled-like chain dynamics, despite the used chain 
length being much shorter than N c . Anomalous scal- 
ing is found even in the short-chain limit. We associate 
the mentioned anomalous features for the Rouse modes 
of the fast component to strong memory effects, which 
are induced by the slow nature of the confining matrix, 
and break the Rouse assumption of time uncorrelation 
of the external forces acting on the tagged chain. 

The model introduces a mixture A/B of bead-spring 
chains with the same number of monomers, N. The 
monomer mass is m = 1. Monomers within a same 
chain are identical (of the same species A or B). We 
simulate systems with N — 4, 7, 10, 15, and 21. The 
monomer- monomer interaction potential is V a p{r) = 
4e[(cr Q/ 3/r) 12 -7c- 12 + 6c" 14 (r/cr Q/3 ) 2 ], where e = 1, c = 
1.15 and a, f3 € {A, B}. Potential and forces are con- 
tinuous at the cutoff r c = ca a p. The interaction diam- 
eters are ctbb = 1, caa = 1-6<7bb, and oab = 1.3<tbb- 
Chain connectivity is introduced by a FENE potential 
EH, V£f NE (r) = -kR 2 eln[l-(Roa aa )-*r 2 ], between 
consecutive monomers, with k = 15 and Rq = 1.5. 
The blend composition is xb — Ab/(Aa + Ab), with 
N a the number of a-monomers. The total number of 
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FIG. 1: Panel (a): Symbols are correlators of the Rouse 
modes for the B-chains, for N — 10 at T — 0.40. Lines are 
fits to a KWW function. Panels (b) and (c): For respec- 
tively the A- and B-chains, p-dependence of the stretching 
exponents, for N = 10 at different temperatures. Symbol 
codes in (b) and (c) are identical. 



pth-mode is g p (i) = f,(t) cas\jjm/N\. Integra- 

tion of the equations of motion yields (X p (t) ■ X g (0)) = 

{2NQ)~ 2 dt'e-^')/^ *V"/M gp (t')g g (*")) 
for the correlators of the Rouse modes. The relaxation 
times are given by t p = (6 2 (/12fc B T) sin -2 [p7r/2iV]. 
At this point no assumption has been made about 
the nature of the forces fj. The Rouse model fully 
neglects spatial and time correlation of the stochastic 
forces, i.e., (g P (t')g q {t")) = l2N(k B T5 pq 6(t' - t"). 
The former two approximations yield respec- 
tively orthogonality and exponential relaxation of 
the Rouse correlators, leading to the expression 
(Xp(t) ■ X,(0)> = (b 2 /8N) sm- 2 \pn/2N]6 pq eM-t/T P ]. 
Orthogonality and exponentiality of the Rouse modes 
are the two main predictions of the Rouse model, 
and are the basis for the derivation of the correlators 
probing chain relaxation [5|. By using the mentioned 
approximations, one obtains Dcu — (NQ~ 1 k B T for 
the diffusivity of the chain center-of-mass [5(. By 
introducing this result in r p (see above) we find the 
scaling relation for the Rouse relaxation times: 



12NDcmt p = b 2 sm- 2 [pir/2N] 



(1) 



monomers is typically 2000-3000. We use x B — 0.3 and 
a packing fraction tf> = 0.53 3]. Temperature T, dis- 
tance, and time t are measured, respectively, in units 
of e/ks, 0BB) and ctbb(«V e) 1 / 2 . Equilibrium simula- 
tions are performed in the microcanonical ensemble, 
with 20-40 independent runs for statistical averages. 
The longest production runs extend to t ~ 4 x 10 6 , 
corresponding to about 800 million time steps. 

Monomer size disparity induces strong dynamic 
asymmetry between the two components by decreas- 
ing T 0. We define the 'dynamic asymmetry', \ = 
TAa/tbBj with r aa the relaxation time of the normal- 
ized density-density correlator F aa (q,t) for a-a pairs 
evaluated at the wavevector q for the maximum of 
the static structure factor Thus, r aa is in the time 
scale of the structural a-relaxation. By using the defi- 
nition F aa (q,T aa ) — 0.3, x increases, from about 4 at 
T = 1.5, more than three decades by decreasing T. 

Now we summarize the main predictions of the 
Rouse model [f|. A tagged chain is represented 
as a string of beads of equal mass connected by 
harmonic springs of constant 3kQT/b 2 , with b the 
bond length. The effective interaction experienced 
by the chain is given by a friction coefficient ( 
and a set of stochastic forces fj. Excluded vol- 
ume interactions are neglected. The chain normal 
(Rouse) modes of index p = 1, ...N — 1 are defined as 
Xp(i) = N^Y^iTjWcoBljpir/N], and follow the 
equations of motion 2A^X p = — fc p X p + g p , with k p = 
2ANk B Tb- 2 sin 2 [ P Tr/2N}. The external force for the 



In the following we test the former predictions in the 
bead-spring blend here investigated. We compute nor- 
malized Rouse correlators for the a-chains, $ pq (t) — 
(X£(t) • X^(0))/([A p 1 (0)] 2 ). According to orthogonal- 
ity 3>p 9 (0) = S pq . For both the A- and B-chains simula- 
tion data fulfill |$p 9 (0)| < 0.1 for all the cases p ^ q at 
all the temperatures, indicating than deviations from 
orthogonality are small. Thus, we conclude that there 
are only weak spatial correlations between the external 
stochastic forces for the chains of both components, 
and Rouse modes as defined above are good normal 
modes. This prediction is not affected by blending, 
and in particular by the observed confinement effects 
on the fast component. Now we test the exponentiality 
of the Rouse modes. Fig. QJi shows results for normal- 
ized Rouse correlators of the B-chains for N = 10 at 
T = 0.40. The decay of & PP (t) from the plateau has 
been fitted [lj| to an empirical Kohlrausch- Williams- 
Watts (KWW) function, cx exp[-(t/T^)^]. Figs. 
and[TJ; show, for respectively the A- and B-chains, the 
T- and p-dependence of the stretching exponent (3 p . Ac- 
cording to the Rouse model relaxation of the Rouse 
correlators is purely exponential {fl p = 1). Clear devia- 
tions from this prediction are observed by decreasing T 
(i.e., by increasing x). At a same T, non-exponentiality 
is systematically more pronounced for the B-chains, 
which show /3 p -values of even 0.3 at T = 0.33. Only 
at high temperature & pp (t) is approximately exponen- 
tial (P p > 0.8) for both components. These features 
arc rather different from observations in the homopoly- 
mer case Q , which even at very low T shows exponents 
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FIG. 2: Panels (a) and (b): Test of Rouse scaling for re- 
spectively the A- and B-chains. Identical symbols corre- 
spond to a same T (triangles: T — 0.33; circles: T — 0.4; 
squares: T = 0.5; diamonds: T — 0.6; crosses: T = 0.75; 
stars: T = 1.0; down triangles: T = 1.5). Identical col- 
ors correspond to a same TV (black: TV = 4; red: TV = 7; 
green: TV = 10; blue: TV = 15; magenta: TV = 21). Solid 
lines indicate power-law behavior ~ (N/p) x . Dashed lines 
correspond to Eq. (JXJ) corrected by a shift factor (see expla- 
nation in Q). For clarity of presentation, data in panel (b) 
are rescaled by a T-dependent factor c. The latter drops 
about two decades from c < 1 at T = 1.5. Panel (c): Static 
Rouse modes (scaled by TV) for A- and B-chains. Lines in- 



dicate an effective power-law 
codes are as in (a) and (b). 



(N/p) ' . Color and symbol 



j3 p > 0.8. Hence, blending induces much stronger non- 
exponentiality for Rouse correlators. 

In principle, non-exponentiality may be related to a 
distribution of intrinsically exponential processes orig- 
inating from structural and/or dynamic heterogeneity. 
It is well-known that a KWW function can be for- 




FIG. 3: Mean squared displacement of middle monomers 
of the B-chains, for TV = 15. Lines indicate an effective 
power-law ~ t y . Filled circles mark the Rouse times n. 



mally expressed as a sum of pure exponential func- 
tions weighted by an adequate distribution G, though 
the latter does not necessarily have physical mean- 
ing. Hence, we express the KWW function for the 



JdT p G(T p ) exp[-t/r } 
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pth-mode as exp[— (t/V^ c )' 3p ] = 
where the distributed values r p are the relaxation 
times of the different exponential processes. Now we 
show that the latter is a good approximation for the 
A-chains but unphysical for the B-ones. In other- 
words , the Rouse modes are basically exponential for 
the A-chains and strongly non-exponential for the B- 
ones. Since the relaxation times for the elemen- 
tary processes would be t p = sin~ 2 (p7r/2TV) (see 
above), with £ a p- independent factor, we can write 
exp[-(t/7*)*>] = J^G(C)exp[-^sin 2 (p7r/2TV)], with 
t = t sin 2 (p7r/2TV) the 'conjugated variable' of £. By 



Tp sin (p7r/2TV), we find exp[ 



■)Pp1 - 



defining f 

J d£G(£) exp[— Since the right side of this equation 
is p- independent, consistency requires the scaling con- 
dition r„ oc sin" ~ 2 (j57r/2TV). Fig. [2] shows a test of Eq. 
([1]) — also in the long wavelength limit t p ~ (N/p) 2 
— for both components at the investigated values of TV 
and T. This comparison also provides a test of the men- 
tioned scaling condition ll|. In agreement with obser- 
vations for the non-entangled homopolymer case in sim- 
ilar bead-spring models [7f or in lattice models [9( , data 
for the A-chains exhibit only small deviations, at small 
N/p. Therefore, exponentiality of the Rouse modes for 
the A-chains is actually a good approximation, and the 
observed non-exponentiality of the Rouse correlators 
basically originates from a distribution of elementary 
exponential processes. The apparently stronger distri- 
bution effects (stronger stretching), as compared to the 
homopolymer [7|, might be related to structural and/or 
dynamic heterogeneities induced by blending. 

The B-chains exhibit a clear breakdown of Rouse 
scaling (Fig. [2Jd). This result demonstrates that for 
the fast B-component, though distribution efects are 
(as for the slow A-component) presumably present, the 
Rouse modes are intrinsically, at low T, of strongly non- 
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exponential nature. Eq. JT]) is approached only in the 
high T limit, where dynamic asymmetry (x ~ 4 at 
T = 1.5) and confinement effects vanish. By decreas- 
ing T we find an ultimate regime r p ~ (N/p) x . 
This crossover is concomitant with the onset of non- 
exponentiality of the Rouse modes. It must be stressed 
that this anomalous scaling behavior is not related to 
particular static features of the Rouse modes. Fig. [2fc 
shows, for all the investigated values of T and N, results 
for the static Rouse modes, (LY p (0)] 2 ). The latter dis- 
play scaling behavior ~ (N/p) 2 - 2 , both for the A- and 
B-chains, in full an alog y with observations for the ho- 
mopolymer case 0, Il2j | . Therefore, anomalous scaling 
in Fig. [2Jd is definitively a non-trivial dynamic effect. 

Exponcntiality of the Rouse modes originates from 
the assumption of time uncorrelation of the stochastic 
forces (see above). Results reported here evidence that 
the latter clearly breaks down for the fast component. 
Several theoretical approaches based on projector oper- 
ator techniques incorporate density fluctuations around 
the tagged chain within a memory kernel in a general- 
ized Langevin equation (GLE). Slow relaxation of the 
memory kernel induces strong non-exponentiality of the 
Rouse modes [13j . The Rouse model (pure exponential- 



ity) is recovered as a particular case for which density 
fluctuations around the tagged chain relax in a micro- 
scopic time scale [l3[ • This rough approximation seems 
to work reasonably for simple non-entangled homopoly- 
mers and for the slow component in the blend model 
here investigated, but it is clearly unrealistic for the fast 
component at low T, which consists of chains diluted 
in a much slower host matrix. 

The mentioned GLE methods, through approxima- 
tions of the memory kernel within renormalized Rouse 
models, predict for t p a crossover (N/p) 2 — > (N/p) x ~ 3,5 
in entangled homopolyme rs 1131) - The latter is in agree- 
ment with simulations 0, Il4l. Il5l Hij of homopolymers 
at fixed T by increasing the chain length beyond the 
entanglement value N c ~ 35. A similar crossover is 
observed in the non-entangled blend model here inves- 
tigated by increasing the dynamic asymmetry between 
the two components [l7| ■ A further analogy in the fast 
B-component with predictions of GLE methods for en- 
tangled homopolymers fl3j is the observed anomalous 
scaling ~ ^a~0-3 f or ^he mean-squared displacement in 
the time interval prior to the Rouse time t\. Fig. [3] 
shows for N = 15 a crossover, concomitant with that 
of Fig. [2t>, from Rouse-like scaling (~ t v ~°- 5 ) to an ul- 
timate anomalous regime as the former. Similar results 
are obtained for the B-chains at the other investigated 
chain lengths. On the contrary, a T-independent expo- 
nent y s» 0.6 is obtained for the slow A-component [3j, 
again in full analogy with the homopolymer case 0] • 

Within GLE methods observation of the mentioned 
anomalous power laws is directly connected to slow 



relaxation of density fluctuations around the tagged 
chain 13]. The latter may be induced by entangle- 
ment, but data reported here for the fast component 
suggest that this is not a necessary ingredient. Analo- 
gies with entangled- like dynamics are indeed observed 
even for N = 4, provided that dynamic asymmetry in 
the blend is sufficiently strong. In summary, results 
presented here suggest a more general frame for chain 
relaxation features usually associated to entanglement 
effects. They also open new possibilities for the appli- 
cation of GLE methods in complex polymer mixtures. 
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